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We report on the first next-to-leading order QCD computation of W + 3-jet production in hadronic 
collisions including all partonic subprocesses. We compare the results with CDF data from the 
Tevatron, and find excellent agreement. The renormalization and factorization scale dependence 
is reduced substantially compared to leading-order calculations. The required one-loop matrix ele- 
ments are computed using on-shell methods, implemented in a numerical program, BlackHat. We 
use the SHERPA package to generate the real-emission contributions and to integrate the various 
contributions over phase space. We use a leading-color (large- A^c) approximation for the virtual part, 
which we confirm in 1, 2-jet production to be valid to within three percent. The present calcu- 
lation demonstrates the utility of on-shell methods for computing next-to-leading-order corrections 
to processes important to physics analyses at the Large Hadron Collider. 
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PACS numbers: 12.38.-t, 12.38.Bx, 13.87.-a, 14.70.-e, 14.70.Fm, 11.15.-q, ll.15.Bt, 11.55.-m 



Particle physicists have long anticipated the discovery 
of new physics beyond the Standard Model at the Large 
Hadron Collider (LHC) at CERN. In many channels, 
discovering, understanding, and measuring new physics 
signals will require quantitatively reliable predictions for 
Standard Model background processes. Next-to- leading 
order (NLO) calculations in perturbative QCD are cru- 
cial to providing such predictions. Leading-order (LO) 
cross sections suffer from large normalization uncertain- 
ties, up to a factor of two in complex processes. NLO cor- 
rections typically reduce the uncertainties to 10-20% l). 

The production of a vector boson in association with 
multiple jets of hadrons is an important process. It forms 
a background to Standard Model processes such as top 
quark production, as well as to searches for supersym- 
metry. Here we present the first NLO computation of 
W + 3-jet production that can be compared directly to 
data, namely CDF results from the Tevatron. 

The development of methods for computing high- 
multiplicity processes at NLO has involved a dedicated 
effort over many years, summarized in ref. [l|. The 
longstanding bottleneck to NLO computations with four 
or more final-state objects — including jets — has been 
in evaluating one-loop (virtual) corrections. Feynman- 
diagram techniques suffer from a rapid growth in com- 
plexity as the number of legs increases; in QCD, NLO 
corrections to processes with four final state objects have 
been limited to the case of all external quarks 0]. On- 
shell methods [1, 0, i, 0, B B E, El [11 , in contrast, do 
not use Feynman diagrams, but rely on the analyticity 
and unitarity of scattering amplitudes to generate new 
amplitudes from previously-computed ones. Such meth- 
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In an on-shell approach, terms in a one-loop ampli- 
tude containing branch cuts are computed by matching 
the unitarity cuts (products of tree amplitudes) with an 
expansion of the amplitude in terms of a basis of scalar 
integrals Q . Recent refinements @, [l^, [ll|, [1^1 , exploit- 
ing complexified loop momenta, greatly enhance the ef- 
fectiveness of generalized (multiple) cuts Evaluating 
the cuts in four dimensions allows the use of compact 
forms for the tree amplitudes which enter as ingredients. 
This procedure drops rational terms, which could be com- 
puted by evaluating the cuts in D dimensions [lij . One 
may also obtain the rational terms using on-shell recur- 
sion, developed by Britto, Cachazo, Feng and Witten at 
tree level [3], and extended to loop level in refs. [1, B|. 

Within the BlackHat program [11 , we determine co- 
efficients of scalar integrals using Forde's analytic ap- 
proach also incorporating elements from the ap- 
proach of Ossola, Papadopoulos and Pittau (OPP) [10|. 
For the rational terms, we have implemented both loop- 
level on-shell recursion and a massive continuation ap- 
proach (related to ZJ-dimensional unitarity) along the 
lines of Badger's method The on-shell recursion 
code is faster at present, so we use it here. The requisite 
speed and numerical stability of BlackHat have been 
validated for one- loop six-, seven- and eight-gluon am- 
plitudes fisll , and for leading-color amplitudes for a vec- 
tor boson with up to five partons 17[, required for the 
present study. (A subsequent computation of one-loop 
matrix elements needed for W ~\- 3-jet production using 
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FIG. 1: Sample diagrams for the seven-point amplitudes qg — + 
evq'gg and qQ — > euq'gQ. The eu pair couples to the quarks 
via a W boson. 



-D-dimensional generalized unitarity within the OPP for- 
malism was described in ref. [l8|.) Other numerical pro- 
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grams along similar lines are presented in refs. 

To speed up the evaluation of the virtual cross sec- 
tion, we make use of a leading-color (large- A^^) approxi- 
mation for the finite parts of the one-loop amplitudes, 
keeping the exact color dependence in all other parts 
of the calculation. Such approximations have long been 
known to be excellent for the four-jet rate in e"'"e~ an- 
nihilation (20| . A similar approximation was used re- 
cently for an investigation of -I- 3-jet production 2ll |. 
which, however, also omitted many partonic subpro- 
cesses. Our study retains all subprocesses. In addition, 
we keep all subleading-color terms in the real-emission 
contributions. In the finite virtual terms of each sub- 
process we drop certain subleading-color contributions. 
"Finite" refers to the e*^ term in the Laurent expan- 
sion of the infrared-divergent one-loop amplitudes in 
e = (4 — D)/2, after extracting a multiplicative factor of 
cr(e) =^(l-He)^2(l-e)/^(l-2e)/(47r)2-^ "Subleading- 
color" refers to the part of the ratio of the virtual terms to 
tree cross section that is suppressed by at least one power 
of either \/N'^ or nj/N^ (virtual quark loops). We mul- 
tiply the surviving, leading-color terms in this ratio back 
by the tree cross section, with its full color dependence. 

For this approximation, we need only the color-ordered 
(primitive) amplitudes in which the W boson is adja- 
cent to the two external quarks forming the quark line 
to which it attaches. Representative Feynman diagrams 
for these primitive amplitudes are shown in fig. [1] Other 
primitive amplitudes have external gluons (or a gluon 
splitting to a QQ pair) attached between the W bo- 
son and the two above-mentioned external quarks; they 
only contribute 22] to the subleading-color terms that 
we drop. As discussed below, we have confirmed that 
for W -\- 1,2-jet production this leading-color approxi- 
mation is valid to within three percent, so we expect cor- 
rections to the W + 3-jet cross-sections from subleading- 
color terms also to be small. 

In addition to the virtual corrections to the cross sec- 
tion provided by BlackHat, the NLO result also re- 
quires the real-emission corrections to the LO process. 
The latter arise from tree-level amplitudes with one ad- 
ditional parton, either an additional gluon, or a quark- 
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TABLE I: Total cross sections in pb for W + n jets with 
^nth-jot ^ 25 GeV as measured by CDF 2]. The resuhs are 
compared to NLO QCD. For W ± 1 and W ± 2 jets, the 
difference between the leading-color approximation and the 
complete NLO result is under three percent. For VK ± 3 jets 
only the LC NLO result is currently available, but we expect a 
similarly small deviation for the full NLO result. Experimen- 
tal statistical, systematic and luminosity uncertainties have 
been combined for the CDF results. 



antiquark pair replacing a gluon. Infrared singularities 
develop when the extra parton momentum is integrated 
over unresolved phase-space regions. They cancel against 
singular terms in the virtual corrections, and against 
counterterms associated with the evolution of parton dis- 
tributions. We use the program AMEGIC+-H ^ to 
implement these cancellations via the C at ani- Seymour 
dipole subtraction method [H]. The SHERPA frame- 
work fisj incorporates AMEGICH — h, making it easy to 
analyze the results and construct a wide variety of dis- 
tributions. For other automated implementations of the 
dipole subtraction method, see refs. (26j . 

The CDF analysis [2,] employs the J ETCLU cone a l- 
gorithm with a cone radius R ~ a/ (Acj))^ + (Ary)^ = 
0.4. However, this algorithm is not generally infrared 
safe at NLO, so we instead use the seedless cone algo- 
rithm SISCONE [28]. In general, at the partonic level we 
expect similar results from any infrared-safe cone algo- 
rithm. For ± 1, 2 jets we have confirmed that distri- 
butions using SISCONE are within a few percent of those 
obtained with the midpoint cone algorithm [29| . 

Both electron and positron final states are counted. 



> 20 GeV, 



> 20 GeV, and Ei^^ > 20 



and the following cuts are imposed 
\tj^\ < LI, $T > 30 GeV, A 
GeV. Here Et is the transverse energy, is the miss- 
ing transverse energy, MjT the transverse mass of the ev 
pair and rj the pseudorapidity. Jets are ordered by Et, 
and are required to have \r]\ < 2. Total cross sections 
are quoted with a tighter jet cut, i?;^* > 25 GeV. CDF 
also imposes a minimum AR between the charged de- 
cay lepton and any jet; the effect of this cut, however, is 
removed by the acceptance corrections. 

CDF compared 0] their measured W + n-jet cross sec- 
tions to LO (matched to partons showers [30|) and the 
then-available NLO theoretical predictions. The LO cal- 
culations differ substantially from the data, especially at 
lower Et, and have large scale-dependence bands. In 
contrast, the NLO calculations for n < 2 jets (using the 
MCFM code 3l| , with the V + 4-parton one-loop matrix 
elements of ref. Q) show much better agreement, and 
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FIG. 2: The measured cross section da{W eu+ > 7i-jets)/di?y* compared to NLO predictions for n = 2,3. In the upper 
panels the NLO distribution is the solid (black) histogram, and CDF data points are the (red) points, whose inner and outer 
error bars denote the statistical and total uncertainties on the measurements. The LO predictions are shown as dashed (blue) 
lines. The lower panels show the distribution normalized to an NLO prediction, the full one for n = 2 and the leading-color 
one for n = 3, in the experimental bins (that is, averaging over several bins in the upper panel). The scale uncertainty bands 
are shaded (gray) for NLO and cross-hatched (brown) for LO. In the n — 2 case, the dotted (black) line shows the ratio of the 
leading-color approximation to the full-color calculation. 



narrow scale-dependence bands. See ref. [2] for details. 

Our aim in this Letter is to extend this comparison 
to n = 3 jets. We apply the same lepton and jet cuts 
as CDF, replacing the cut by one on the neutrino 
Et, and ignoring the lepton-jet AR cut removed by 
acceptance. We approximate the Cabibbo-Kobayashi- 
Maskawa matrix by the unit matrix, express the W cou- 
pling to fermions using the Standard Model parame- 
ters OfQED — 1/128.802 and sin^ 9w — 0.230, and use 
mw = 80.419 GeV and Tw = 2.06 GeV. We use the 
CTEQ6M [ail parton distribution functions (PDFs) and 
an event-by-event common renormalization and factor- 
ization scale, fi = -|- p|,(VF). To estimate the scale 
dependence we choose five values in the range (^, 2) x /z. 
The numerical integration errors are on the order of a 
half percent. We do not include PDF uncertainties. For 
T4^ + 1, 2-jet production these uncertainties have been 
estimated in ref. In general they are smaller than 
the scale uncertainties at low Et but larger at high Et- 
The LO calculation uses the CTEQ6L1 PDF set. For 
n = 1,2 jets, NLO total cross sections agree with those 
from MCFM 31|, for various cuts. As our calculation is 
a parton-level one, we do not apply corrections due to 
non-perturbative effects such as induced by the underly- 
ing event or hadronization. Such corrections are expected 
to be under ten percent 2]. 

In table U we collect the results for the total cross 
section, comparing CDF data to the NLO theoretical 



predictions computed using BlackHat and SHERPA. 
The columns labeled "LC NLO" and "NLO" show respec- 
tively the results for our leading-color approximation to 
NLO, and for the full NLO calculation. The leading-color 
NLO and full NLO cross-sections for W ^ \- and 14^-1-2- 
jet production agree to within three percent. We thus 
expect only a small change in the results for W -I- 3-jet 
production once the missing subleading-color contribu- 
tions are incorporated. 

We have also compared the Et distribution of the n*'' 
jet in CDF data to the NLO predictions for W + 1,2, 3- 
jet production. For W + 2, 3-jets these comparisons are 
shown in fig. [2l including scale-dependence bands ob- 
tained as described above. For reference, we also show 
the LO distributions and corresponding scale-dependence 
band. (The calculations matching to parton showers (30j 
used in ref. Q make different choices for the scale varia- 
tion and are not directly comparable to the parton-level 
predictions shown here.) The NLO predictions match 
the data very well, and uniformly in all but the high- 
est Et bin. The central values of the LO predictions, 
in contrast, have different shapes from the data. The 
scale dependence of the NLO predictions is substantially 
smaller than that of the LO ones. In the W + 2-jet case, 
we also show the ratio of the leading-color approxima- 
tion to the full-color result within the NLO calculation: 
the two results differ by less than three percent over the 
entire transverse energy range, considerably smaller than 
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FIG. 3: The theoretical prediction for the Ht distribution in 
W + 3-jet production. The curves and bands are labeled as 
in fig. H 



the scale dependence (and experimental uncertainties). 

In fig. [21 we show the distribution for the total trans- 
verse energy Ht, given by the scalar sum of the jet and 
lepton transverse energies, Rt = Ylij ^T^j + + ^t- 
We show the NLO and LO predictions, along with their 
scale-uncertainty bands. As in the distributions, the 
NLO band is much narrower; and the shape of the dis- 
tribution is altered at NLO from the LO prediction. 

In summary, we have presented the first phenomeno- 
logically useful NLO study of -I- 3-jet production, and 
compared the total cross section and the jet Et distribu- 
tion to Tevatron data The results demonstrate the 
utility of the on-shell method and its numerical imple- 
mentation in the BlackHat code for NLO computations 
of phcnomcnologically- important processes at the LHC. 
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